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Stochastic  Duels--II 
by 

Trevor  Williams 

System  Development  Corporation 
Santa  Monica^  California 

ABSTRACT 

This  report  continues  the  project  of  systematically  developing  the  elements 
of  a  stochastic  theory  of  combat  begun  in  an  earlier  work.  The  mean^  variance^ 
third  central  moment^  etc.,  of  a  marksman’s  time  to  hit  a  passive  target  are 
evaluated  in  terms  of  the  corresponding  parameters  of  his  time  to  fire  a  single 
round.  Next,  the  solution  of  the  simple  duel,  in  the  case  where  each  protago¬ 
nist’s  time-to-kill  is  distributed  as  a  gamma-variate,  is  obtained  as  the  cumu¬ 
lative  distribution  of  a  certain  binomial  variate,  and  this  result  is  employed 
to  furnish  an  approximate  solution  to  the  general  simple  duel.  An  expansion  of 
the  moment -genera ting  function  of  the  marksman's  time-to-kill  in  powers  of  his 
kill  probability  is  next  derived  and  found  to  provide  a  good  approximation  to 
the  solution  of  the  simple  duel;  various  properties  of  the  expansion  are  also 
considered.  We  then  examine  a  stochastic  battle  in  which  all  men  on  both  sides 
are  at  all  times  able  to  participate  in  the  action;  it  is  deduced  that  for  large 
numbers  of  opposing  forces  the  lead-order  behavior  of  the  system  is  deterministic 
and  reduces  to  Lanchester’s  square  law,  the  exchange  rate  being  suitably  interpreted. 
Finally,  we  entertain  a  stochastic  battle  where  the  two  forces  can  only  be  brought 
into  play  at  a  single  point  of  contact:  we  find  that  the  lead -order  behavior  is 
now  given  by  Lanchester’s  linear  law,  but  in  this  case  we  are  able  to  go  much 
further  toward  obtaining  a  complete  solution. 
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1.  INTRODUCTION 

After  a  tiresom  harangue  in  his  ordinary  style^  he  took  down 
from  his  book  shelves  a  number  of  musty  volumes  on  the  subject 
of  the  duello,  and  entertained  me  for  a  long  time  with  their 
contents.  .  . 

---  Poe  (Mystif ication) 

The  present  report  is  a  continuation  and  extension  of  the  project,  begun 
by  Williams  and  Ancker  in  an  earlier  paper,  of  developing  the  elements  of  a 
stochastic  theory  of  combat  by  means  of  idealized  mathematical  models  whose 
implications  can  be  thoroughly  traced  out.  Free  use  is  made  of  the  results  of 
that  prior  work,  which  we  henceforth  refer  to  as  Stochastic  Ducls--I:  an 
equation  cited  here  as,  say,  "(1-22)"  means  the  22nd  equation  of  that  paper, 
and  the  citation  "I-Figure  3"  means  the  third  figure  of  that  paper.  The  term 
"stochastic  duel"  is  a  convenient  one  and  we  now  broaden  it  to  include  any 
situation  in  which  two  enemy  forces  engage  in  conflict;  in  this  we  have  a  certain 
etymological  authority  since  the  Latin  "bellum"  is  itself  a  variant  form  of  the 
older  "duellum, "  which  meant  a  contest  between  two  parties,  and  gives  us  our 
own  "duel"  by  way  of  the  Italian  "duello. "  The  discipline  of  model-building 
in  any  application  of  mathematics  to  the  physical  world  requires  certain 
qualifications  if  it  is  to  be  successful.  First  of  all,  significant  aspects  of 
the  problem  must  be  ferreted  out  for  consideration  and  in  so  doing  we  must  hew 
away  irrelevancies.  There  are  many  instances  in  the  literature  of  military 
operations  research  where  the  argument  is  cluttered  up  with  minute  details  which 
lend  a  spurious  verisimilitude  to  the  model  but  have  no  bearing  on  the  fundamental 
issues  and  render  the  mathematics  completely  intractable.  Having  sorted  out 
the  key  features  of  the  system,  we  must  weigh  what  we  might  like  to  do  against 
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what  we  can  reasonably  expect  to  be  able  to  do.  An  idealized  and  simplified 
model  which  we  can  solve  explicitly  and  completely  is  of  much  greater  utility 
in  fostering  insight  than  a  very  realistic  model  which  we  cannot  solve.  Finally^ 
we  must  build  up  the  bits  and  pieces  of  information  so  obtained  into  more  and 
more  comprehensive  patterns  in  order  to  approach  the  complexity  of  the  world  we 
live  in.  This  report  represents  a  modest  attempt  to  take  a  fev7  such  steps. 


2.  CONNECTION  BETWEEN  CUMULANTS  OF  TIME-TO-FIRE  AND  TIME-TO-KILL 

And  a  certain  man  drew  a  bow  at  a  venture^  and  smote  the  king 
of  Israel  between  the  joints  of  the  harness.  .  . 

- 1  Kings  22;  34 


We  assume  that  a  marksman  fires  repeatedly  at  a  target  until  he  hits  it^ 
his  probability  of  doing  so  on  any  given  round  being  p,  where  p  is  constant 
from  round  to  round.  The  frequency  function  of  his  time  to  fire  we  call  f(t) 
and  that  of  his  time  to  score  a  hit,  h(t).  The  corresponding  moment -genera ting 
functions  we  write  as 


9(z) 


w 

-L 


e^^f (t)dt 


and  il>(z) 


e^*'h(t)dt 


(2.1) 


Assuming  successive  firing-times  to  be  statistically  independent,  we  have  shown 
(1-14)  that 


pa>(2) 

l-qcp(z) 


(2.2) 
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Strictly  speaking,  this  result  was  derived  for  the  characteristic  function 
rather  than  the  moment-generating  function,  but  since  one  may  be  obtained  from 
the  other  by  the  subsituation  z=iu,  i. e. ,  a  rotation  through  90°,  the  distinction 
is  of  no  importance  analytically. 

The  cumulant-generating  function  for  the  time-to-fire  will  be  called 

00 

e(z)  =  log  cp(z)  =  Z)  nl  ’  (2-3) 

n=l 


where  is  the  mean  time-to-fire,  the  variance,  etc.  Similarly,  for  the 
cumulant-generating  function  of  the  time-to-kill  we  write 


00  n 

©(z)  =  log  <b{z)  =  Z)  K  ^ 


n=l 


n  nl 


(2.4) 


Then  (2. 2)  yields 


-qcp 


-qe 


(2.5) 


upon  invoking  (2.3);  in  other  words. 


1  ®  o 

0  =  0.log  =  e-iog[i  -  ^  (e^-i)] 


(2.6) 


For  convenience  writing 


t  =  q/p  ^ 


(2.  7) 


(2.  6)  becomes 
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00  V  “  V  V 

G  -  0  +  S  f  (e®-!)''  -  0  +  S  ^ 

v=l  v=l  ^=o 


00  V  V  00  n^n 


V  =  1  (1=0 


00  .n  00  .V  V 


e+Efr  Ef 

n=o  v=l  p=o 


(2.8) 


But^  we  recall  that^  by  definition 


=  \+r“k 


(2.9) 


o‘'\  -  oca'-^Ui^) 


(v=2^  3_, . .  .  ) 


(2.  10) 


from  which  it  follows  that 


A  ^  X"'  /•  \  ”  “h  A  ►  \ 


(2.11) 


(See^  e.  g,  ^  Whittaker  and  Robinson.)  Hence^  for  the  sequence  Uj^=k  we  have 


1=0  J  k=o 


(2.12) 


the  latter  form  being  the  conventional  abbreviation.  Since  differences  of  a 
polynomial  of  higher  order  than  its  degree  vanish,  we  have 
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=  0 


(v  >  n) 


(2.  13) 


and  substituting  this  result  and  (2.11)  into  (2.8)^  we  obtain 


“  ^n  11  y 

©  =  6  +  E  ^  E  —  . 

^  nl  1  V 
n=o  v=l 


(2.  14) 


We  note  that  the  term  for  n=o  is  a  null  sum  and  may  therefore  be  suppressed. 
Furthermore^  we  have 


fiV  -  v(aV^  +A''-V-b 


(2.  15) 


For^  in  consequence  of  (2.9)  and  (2.10)^ 


(2.  16) 


hence  the  right-hand  side  of  (2. 15)  is 


1  1  v-1  T  T  1 

.v-l,n-l  .v-l-n,v-l.  .  .inH-I 

V  .  A  1  =  V  E(-)  (  H  )  (4+1)  7 

^L=0 


(2.  17) 


by  (2.11)^  and  this  last  may  be  re-written  as 


V 

V*/  ^v-M.,v-l.  n-1  'o,  sV-u  V  /V-1.  n 

V  Z7(-)  "^(^-1)4  =  2_/(-)  ~  (n-l)4 


4=1 


4=1 


=  E(-)"-"(^)4"  =  aV  , 

4=1 


(2.  18) 


when  we  refer  back  to  (2.  12)  and  note  that  the  term  for  n=o  contributes  nothing 


13  September  1963 


8 


SP-1017/003/00 


to  the  sum.  Hence  (2.  15)  is  established;  when  we  employ  it  in  (2. 14)  we  obtain 


0=0+  S 

n=l 


£1 

nl 


v-l 


+  a'-V’S 


00 

=  0  +  s 

n=l 


n 


n 


‘V=l 


+'St''+Vo"-i 

V=0 


(2.  19) 


Now^  of  the  sums  within  the  parentheses^  the  latter  may  be  extended  up  to  v=n^ 

since  this  entails  merely  the  introduction  of  a  term  t'^A’^o'^  ^  which  vanishes 

by  virtue  of  (2.13);  also  the  former  may  be  extended  down  to  v=o^  introducing  a 
n  1 

term  0  j  which  also  vanishes  provided  n  >  1.  Hence  we  may  rewrite  (2. 19)  as 

e  =  0  +  at  +  i:  (L+t)  EtW-‘ 

n=2  ^  v=o 


n=l 


n 


=  (1+t)  E  ^ 


v=o 


(2.  20) 


upon  noting  that  the  term  for  n=l  precisely  compensates  for  the  term  in  front 
of  the  summation.  From  (2.7)^ 

1+t  =  1/p  ,  (2.21) 


and  thus^  writing  (2.20)  out  at  some  length,  in  terms  of  the  well-known  values 
of  the  "differences  of  zero,"  we  obtain,  finally. 


0  =  70+\'^+'^  (1+2  5)17  +  ^(1+65  +  +  •  •  • 


2  21  '  2  p'3I  '  „2^‘'''p  '  “„2 

p  p  P  P 


(2.  22) 
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Going  back  now  to  (2.3)  and  (2.4)  and  substituting  these  expansions  into 
(2.22)  and  comparing  coefficients  of  powers  of  z,  we  readily  derive  the  following 
relationships  between  the  K’s^  which  are  the  cumulants  of  the  time-to-kill^  and 
the  K'S;  which  are  the  cumulants  of  the  time-to-f ire: 


and 


K. 


K, 


1  /  ,2. 

=  ^[p^k2+3pq/c^/c2+q(i+q)'^ij 


> 


(2.  23) 

(2.24) 

(2.25) 


=  ^[p\^+p^q(4kj^/C2+3K2)+6pq(l+q)KlK2+q(p^-l-6pq+6q^)/C^]  .  (2.26) 


These  are  the  desired  results.  A  convenient  check  on  the  algebra  is  afforded 
by  the  fact  that  for  exponentially-distributed  firing  time  with  rate  of  fire  r, 
we  have 


cp(z) 


00 


1 

1-z/r 


(2.27) 


by  (2.  1)  and 


'I>(z) 


2 _  =  1 

z  1-z/rp 

7  ■  'I 


(2.  28) 


by  (2.2).  Also^  from  (2.3)  and  (2.27)  it  follows  that 
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n 


-  n  n  1 
11=1 


-  iog(i  -  f)  =  E 


n 


^1 


n=l 


(2. 29) 


whence 


K 

n 


(n-1)  ! 


(2.  30) 


and  similarly^ 


K 

n 


(n-1) I 
n  n 
P 


(2.  31) 


Substitutions  into  (2. 23) - (2. 26)  yields 


and 


i-  =  i  .  i 

rp  p  ’  r 


(2.  32) 


1  _  JL/  ^ 

2  2  ~  2'^  2  ^  2^ 

r  p  p  r  r 


(2.  33) 


^  +  3pq  -^  +  q(l+q)  ~] 
r'^p  p  r  r"^ 


(2.  34) 


6  lr36^2.8^3v^,  ,i^vl 

”7^4  "  “4^P  ~  ^  P  ^  6pq(i+q)  — 

r  p  p  r  r  r  r 


2  2  1 
+q(p  +6pq+6q  )-7'l  . 


(2.  35) 


The  first  two  of  these  are  obvious  identities;  the  latter  pair  reduce  to 
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2p^+3pq+q(l+q)  =  2  , 

and 

3  2  2  2 

6p  +llp  q+6pq(l+q)+q(p  +6pq+6i  )  =  6  ^ 

which  may  be  verified  by  replacing  p  by  1-q  throughout. 


(2.  36) 


(2.  37) 


3.  SOLUTION  OF  THE  SIMPLE  DUEL  WHEN  THE  TIMES -TO-KILL  ARE  GAMMA- VARIATES 

Laertes.  My  lord^  I'll  hit  him  now. 

King.  I  do  not  think 't. 

- Hamlet  (V,  ii) 


The  solution  of  the  simple  duel  which  we  have  already  exhibited  (1-20) 
becomes^  in  terms  of  moment-generating  functions^ 


P(A) 


e+i«> 


e-ioo 


(3.  1) 


As  we  might  expect,  the  only  difference  in  the  result  is  that  the  path  of 

O 

integration  has  been  rotated  through  90  in  the  complex  plane.  Now,  we  define 
a  7(k,a)  variate  as  one  having  as  frequency  function 


f(t) 


,k-l 


(k-1)  la 


-t/a 


(3.2) 


The  MGF  (moment-generating  function)  is  thus 


<I>(z) 


k-1  , 

- r 


(k-1) la 


1 

(l-az)*^ 


(3.  3) 
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This  result  is  obvious  v^en  we  recall  that  a  /(k^a)  variate  is  the  sum  of  k 
independent  exponential  variates^  each  with  mean  OU 

Let  us  now  assume  that  A’s  time  to  kill  is  distributed  as  7(k^0:)^  while 
B's  is  7(£^P).  Then  we  have 


P(A) 


e+ioo 

JL  r  - i - 


(3.  4) 


This  is^  of  course^  merely  Prob(7(k^a)  <7(i^P))-  Writing 


=  a/p  } 


and  replacing  z  by  z/p  in  (3. 4)^  we  find 

€+i<» 

1 


dz 


G-loo 


(l+4-z)^(l-z)^  ^ 


(3.5) 


(3.6) 


Strictly  speaking^  the  e  in  (3.6)  is  not  the  same  as  that  in  (3.4)^ 
need  be  concerned  with  is  the  fact  that  it  is  some  positive  number, 
the  contour  in  (3.6)  into  the  right  half-plane  since  by  doing  so  we 
residue  (namely  the  one  at  z=l)  rather  than  two.  Since  the  contour 
the  pole  in  the  negative  sense^  we  obtain 


but  all  we 
We  complete 
pick  up  one 
encircles 


P(A) 


res  _ i _ 

(l+\lfz)^(l-z)’^z 


(3.  7) 


Writing 
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this  becomes 


P(A)  =  -  coeff  -  1 - - 

^  (l+t+'ifv;)  (-w)  (1+w) 

.  coeff  _ I _ 

(1+\1/)'^  (1+  ■^^)’^(l+w) 


£-1 

E 

v=o 


_ 1 

(l+4r) 


k+i-1 


i-1  ,  .  1 

Zy(-)  (v)'l'  (1+t) 


,  l-l  ,  l-v-l 

- V?n  E  ( 

(l+\}f)  v=o  (1=0 


Z-v 


E(-)^t)  E(V-;V 

(l+i]/)  v=o  |a=v 


. \4, 1  sV  B(-)‘'(;‘')(7-vb 

(1+t)  ^  (1=0  v=o 


(3.9) 


The  inner  summation  here  may  be  evaluated  by  means  of  two  elementary  and 
well-known  identities,  namely. 


*  -  (-)''(’'Tb 


(3.  10) 


and 
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S(5>(t.5:v)  -  <T) 


v=o 


(3.  11) 


The  first  of  these  is  self-evident^  and  the  second,  known  as  the  "factorial 
binomial  theorem,"  and  due  to  the  mathematican  Vandermonde,  may  be  obtained  by 
equating  coefficients  of  like  powers  of  t  in  the  identity 


(3.  12) 


Using  (3.10)  twice  and  (3.11)  once,  the  inner  sum  in  (3.9)  may  be  written 


c-)^  E(?)(t;:v-)  = 


v=o 


(3. 13) 


and  so  we  arrive  at  the  result 


(1+t)  M=o 


(a+p) 


(3.  14) 


upon  invoking  (3.5).  In  other  words,  if  we  set 


(3.  15) 


and  define  to  be  a  binomial  variate. 


[1  B(n,p) 


(3.  16) 


then 
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P(A)  =  Prob(pi  <  i-l)3  (3.  17) 

may  be  looked  up  in  the  standard  tables.  We  may  summarize  the  result  of  our 
calculation  in  the  single  identity 

Prob(7(k,a)  <  r(i.P))  =  Prob(B(k+i-l,  <  £-l]  .  (3.18) 

We  may  combine  this  result  with  those  of  the  first  section  to  derive  an 

approximate  solution  to  the  simple  duel.  It  follows  from  (3.3),  and  is  of  course 

2 

well-known,  that  for  a  7(k,0!)  variate  the  mean  is  ka  and  the  variance  ka  . 
Identifying  these  with  corresponding  cumulants  of  the  time-to-kill  in  (2.4),  we 
have 


=  ka  and 


(3.  19) 


from  which  it  follows  that 


2 

a  =  ^2^^!  ^  ~ 

We  now  make  use  of  (2.23)  and  (2.24)  and  make  the  fact  that  the  quantities  k 
and  a  refer  to  side  A  explicit  by  putting  a  subscript  A  on  all  the  symbols; 

a  =  (Pa'"1+<1a^1)  Va  ^  =  ^^a/^PaVVa)  • 

2 

Here  we  have  also  replaced  and  by  the  commoner  n  and  a  .  Likewise  we 


have  the  pair  of  relations 
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P  =  ^  "  Pb/^Pb'^b'^bpI)  •  <3.22) 

Utilizing  (3.21)  and  (3.22)  in  (3.  15) -(3.  17)  gives  rise  to  the  approximation 
alluded  to.  Of  course,  k  and  i  are  no  longer  integers  here  and  interpolation 
in  the  tables  of  the  binomial  distribution  becomes  necessary. 

Certain  limiting  cases  are  of  interest.  First  of  all,  we  may  dismiss  the 
situation  where  p^-»-  1,  for  then  (3.21)  merely  falls  back  into  (3.20)  and  indeed 
all  the  K's  become  equal  to  the  corresponding  /c’s,  as  they  must,  since  the  time- 
to-kill  becomes  identical  with  the  time-to-f ire.  The  same  is  true  when  p^->  1. 

D 

When,  however,  p^  and  p^-*-  o  we  find  from  (3.21)  and  (3.22)  that 

P'B 

^  ~  k>l,  i^l  (p  ,  p  ■+•  o)  .  (3. 23) 

Pa  Pb  a  ^ 

Now,  in  (I-l),  we  defined 


(3.  24) 


Hence  in  this  limiting  case  we  have,  from  (3.17), 


P(A)  -  Prob(B(l,  <  0)  =  ^  = 


Pa'a 


Pa'a-^Pb'b 


(3.25) 


which  is  the  same  as  (1-6).  Indeed  the  left-hand  side  of  (3.18)  tells  us  that 
we  are  dealing  with  precisely  the  same  problem,  since,  with  k=i“l,  both  gamma - 
variates  are  merely  exponential-variates.  Equation  (3. 25)  suggests  that  we 
might  expect  our  isoprobs  to  radiate  from  the  origin  of  the  (p  , p  )-plane  the 

/V  o 
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same  way  as  in  Figure  2  of  Stochastic  Duels-I^  and  we  shall  subsequently  find 
further  confirmation  of  this. 

Another  case  which  recommends  itself  is  that  where  the  coefficients  of 
variation  of  A's  and  B’s  time-to-fire  are  small;  we  might  at  least  hope  that 
this  would  be  a  reasonable  assumption  in  practice.  We  then  find  that 


a  -o 


k  — 


^B 


where  the  symbol  "»"  carries  the  customary  meaning  of  "is  much  greater  than.  " 

If  we  further  assume  that  q.  and  q  are  small^  as  many  military  men  would  insist^ 
then  (3.16)  approaches  normality  and  (3.17)  may  be  evaluated  accordingly. 

An  illustration  may  help  to  clarify  the  relations  between  the  approximations 
we  have  obtained.  In  (1-24)  we  have  derived  the  solution  to  the  simple  duel 
with  firing-times  distributed  as  gamma-variates  with  k=2.  For  this  case  we  have 


2  12 
^A  ""  2^A 


and 


2  12 

'^B  =Fb 


(3.27) 


Hence  (3.21)  and  (3.22)  become 


a  =  ^  ,  P  = 


2"aPa 


2Vb  ^ 


2  2 
k  =  — -  and  Z  = 


1+9t 


For  simplicity^  take 


=  Pb  =  3  and  =  2tjj  ; 


(3.  28) 


(3.29) 


then  the  exact  solution^  (1-24),  yields 
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P(A)  =  -fl  =  •  71795  ,  (3.  30) 

whereas  (3.17)  leads  to 

P(A)  «  Prob(B(2,^)  <  |)  .  (3.31) 


Now,  when  /  is  an  integer 

,  all  terms  up  to  and  including  i-1  are  to  be  taken  in 

evaluating  (3.17);  hence 

we  may  display  the  complete 

cumulative  probability 

distribution  implicit  in 

(3.31)  as  follows: 

£ 

C.  D.  F  A 

A^ 

-1 

0 

4/9 

0 

4/9  0 

4/9 

-1/3 

1 

8/9  -1/3 

1/9 

2 

1 

Table  1. 


We  have  here  taken  the  range  of  i  to  be  such  that  the  first  column  of  differences 
consists  of  the  entire  binomial  distribution  of  order  2,  with  p  =  1/3;  we  thus 
retain  all  the  information  at  our  disposal.  First  of  all,  linear  interpolation 
in  Table  1  yields 

2 

P(A)  ,  7.14%  smaller  than  true  value.  (3.32) 
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It  is  interesting  to  remark  that  this  is  the  same  answer  we  would  get  upon 
substituting  (3.29)  into  (1-6),  the  solution  of  the  simple  duel  with  exponential 
firing-times.  That  this  might  provide  a  fair  approximation  may  be  seen  from  the 
fact  that  the  isoprobs  of  I-Figure  5  are  all  reasonably  straight.  Applying 
the  Newton-Bessel  interpolation  formula  to  Table  1,  on  the  other  hand,  we  find 

P(A)  «  I  +  ^  ^  '  2  3^  "  i|  "  ^  small.  (3.33) 

Hence  in  this  case  the  use  of  the  more  refined  approximation  has  not  quite  cut 
the  error  in  half,  and  furthermore  remains  on  the  same  side  of  the  true  value. 

It  seems  likely,  however,  that  duels  for  which  the  isoprobs  display  greater 
curvature  would  give  rise  to  a  more  consistently  favorable  behavior  of  (3. 17) 
over  that  of  (1-6). 

4.  APPROXIMATION  TO  THE  SOLUTION  OF  THE  SIMPLE  DUEL  IN 
TERMS  OF  MOMENTS  OF  THE  TIMES -TO -FIRE 


Round  numbers  are  always  false. 


- Dr.  Johnson 


Looking  at  (2. 23) -(2.  26)  we  see  that  if  we  define  a  new  variate 


T  =  pt  , 


where  t  is  the  time-to-kill,  then 


(4.  1) 
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(n-1)  l/<"  (p->o)  ,  (4.  2) 

tVi 

where  K  (t)  is  the  n—  cumulant  of  the  new  variate  t  and  without  any  argument 
n  1 

is  the  mean  of  the  time-to-f ire.  But,  it  follows  from  (4.  2)  that  T  is  ex¬ 
ponentially  distributed,  with  mean  i.  e.  , 

T  'v-  7(1, ,  (4.  3) 

or,  in  virtue  of  (4.1), 

t  ~  7(1^  ~)  •  (^- 

Hence,  when  p^  and  p^  are  both  small,  the  times-to-kill  are,  in  the  limit, 
exponentially  distributed,  with  means  3nd  p^/p^,  respectively.  The 

solution  therefore  reduces  to  the  one  for  exponential  firing-times,  for  which 
we  saw  (1-2)  that  the  situation  just  described  holds  rigorously. 

This,  then,  yields  the  zeroth-order  behavior  of  the  solution  and  raises 
the  natural  question  of  whether  one  can  push  the  analysis  further  to  obtain 
higher-order  corrections,  thereby  taking  into  account  the  curvature  of  the 
isoprobs  for  all  but  the  case  of  exponential  firing-times.  The  multiplication 
of  the  time-to-kill  by  p  in  equation  (4. 1)  results  in  a  random  variable  with 
MGF 


$^(z)  =  E(e’"^)  =  E(eP‘^^)  =  ^^.(pz)  .  (4.5) 

!,  Hence,  by  (2.2),  we  are  led  to  consider 

1 

c 

i 


¥ 
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<t(pz) 


pq>(pz) 

l-qcp(pz) 


(4.6) 


Employing  i’Hopital’s  rule  we  find 


lim 

p-vo 


<I)(pz) 


cp(o) 


lim  _ E _ 

p-*-o  l-(l-p)cp(pz) 


lim 

p->-o 


_ 1 _ 

cp(pz)-(l-p)zcp''(pz) 


1 


1-zcp  (o) 


1-pz 


(4.  7) 


where  is  the  mean  time-to-fire.  We  recognize  (4.  7)  at  once  as  the  MGF  of  an 
exponential  variate  wi''h  mean  and  thus  we  have  a  corroboration  of  (4.3). 
Prompted  by  the  last  result  v'e  now  write 

00  k 

<I>(pz)  =  Z/\(z)  )  (4.8) 

k=o 

and  recall  that 


00  k  k 

<P(p^.)  -  ^  . 

k=o 


(4.  9) 


where  we  use  to  designate  the  k —  moment  about  the  origin^  of  the  time-to- 
fire.  Equation  (4.6)  then  yields 


P 

k=o 


k  k 
kl 


(4.  10) 


The  term  within  the  accolades  is  merely 


00  k  k  00  k+1  k 

o  o 


=  Ep' 
1 


k-1 


kn 


^k-1  (k-1) 1  ■  ^k  kl 


(4.11) 
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and  so  (4. 10)  becomes 


1 

k=l 


v=o  k=l  '  ' 


•1  k 


But  the  left-hand  side  of  this  equation  may  be  rewritten  as 


k= 


v=l  '  ' 


V  V 


v=l  k=l 


(4.  12) 


(4.  13) 


upon  inverting  the  order  of  summation.  Comparing  the  coefficients  of  like  powers 
of  p  in  (4.12)  and  (4.13)^  we  find 

S(k)z*^'^(l9ij^_^-zHj^)\|f^_j^(z)  =  (v=l,2;3;...)  .  (4.14) 

k“  1 

Repeated  application  of  this  formula  enables  us  to  solve  successively  for  the 
\|f's.  We  readily  find 

=  1^  ^  15) 

(p  -2p^)z^ 

t,(z)=— - —  ,  (4.16) 

2(l-^z)^ 

t,(z)  =  - - - j  {(^„-6^n  +6^5)z^  +  y(3|i^-2nn  )z'^)  ,  (4.17) 

^  3(l-pz)^  ^  ^23 


and 
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'1^3  (z)  =  - - - 3 

4(1-m.z) 

2  2  5 

+(8m.  m.2-I-4|12M-3-12m.m.2-2M-M.4)z 

+(3|i2-W2M-3+!^^M-4)z^}  •  (^-18) 

From  here  on  the  complications  become  well-nigh  insuperable.  It  is  very 
questionable^  however^  whether  we  would  ever  want  to  go  any  further  than 
since  \|r^  introduces  the  fifth  moment  of  the  time-to-f ire;  all  we  ever  really 
know  about  these  moments  are  statistical  estimates  and  we  recall  that  the  sampling 
variance  of  the  k~  moment  entails  moments  up  to  and  including  the  (2k)—. 

Hence^  it  is  extremely  difficult  to  get  any  kind  of  hold  on  the  higher-order 
moments  of  a  distribution  unless  the  sample  size  is  very  large  indeed. 

A  check  on  the  algebra  of  the  last  four  equations  is  afforded  by  considering 
once  again  the  case  of  exponential  firing-times.  There  we  have 

=  nlp'^  (n=0, 1,2,.  .  .  )  ,  (4.19) 


and  when  we  substitute  these  values  into  (4.  15) -(4.  18) ,  we  see  that  >^2  and 
vanish  identically.  This,  of  course,  re-af firms  equation  (1-2). 

We  shall  have  a  little  more  to  say  about  the  \|> '  s  later  in  this  section,  but 
for  the  present  it  will  suffice  to  note  that  their  general  form  is  obviously 


♦„(»)  - 


~TT  Zyt  z 

n+1  nv 


n+v 


(l-pz)“'"  v=o 


(4.  20) 


(The  doubly-subscr ipted  will  not  be  confused  with  the  singly-subscripted  and  the 


13  September  1963 


24 


SP-1017/003/00 


dual  use  of  presents  a  needless  proliferation  of  the  notation.)  We  note  that 


\1/  is  iust  a  Kronecker  delta: 

^no 

=  l(n=o)  ^  =  0  (otherwise).  (4.21) 

The  following  \^'s  are,  by  (4.  16) -(4. 18) , 

til  = 

^21  =  ■|(fi3-6pp2+6fi^)  ,  (4.23) 

^22  =  i<3P2-2fJi^3)  ,  (4.24) 

^31  =  •|(H4-8ii|i3+36p^|i2-6p2-2A|i^)  ,  (^-25) 

2  2  1 

^22  =  2p  P2+^l2^^3■3l^^^2  ‘  2^^4  ^  (^*26) 

and 

'*^33  "  i^^2'‘^^2^3'^^^4^  •  (^-27) 
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and 


1 


k+v 


,k+l  ^"kv 
(l-pigZ)  v=o 


(4.31) 


Here^  e. g. ^  we  clearly  have 


♦llW  -  < 


(4.32) 


where  is  the  second  moment  about  the  origin  of  A's  time-to-f ire;  and  so  on. 

^  A 

Substituting  into  (3.1)  we  obtain,  upon  invoking  (4.28)  and  (4.29), 
e+i"  00 


k 

p  oo 


/p  W  p 

f;;)  kf  ^ 

e-ioo 


£ 

z 


(4.  33) 


Writing  for  convenience 


a  =  p^/p^  and  p  =  pg/pg  , 


(4.  34) 


(4. 30)  and  (4. 31)  now  yield 

e+i® 


P(A)  = 


1 

2rti 


k! 


£ 

oo  n 

_  V  ^ 

.k+1  “  £' 


e-ioo 


k=o  (l-Kfz)''""  £=o  (1-pz) 


£+1 


v=o 


X=o 


iX' 


(4.  35) 


If  now  we  write 


^Ki; 

"vX 


e+i® 

r 

2jti./ 

e-ioo 


k+;!+v4-X-l  , 
z _ dz 

,1  .k+l,,  ^  .  ^+1  ■’ 

(l+f/z)  (1-Pz) 


(4.  36) 
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then  (4. 35)  becomes 


P(A)  -  S  ^  S 

k=o  i=o  v=o 


k+v 


\,(A) 


i 

E 

X=o 


Vi 


(4.37) 


Now^  when  we  attempt  to  evaluate  this  we  are  immediately  confronted  by  the 
unpleasant  fact  that  the  integral  in  (4.36)  only  converges  when  v  +  X  <  1. 

It  is  easy  to  calculate  by  residue  theory  that 


=  -g- 

oo  QH-p  ^  10 


1 

2  ^ 

(.c^y 


and 


1 

2  ^ 

(a4^) 


(4. 38) 


but  as  soon  as  we  consider  the  next  case,  where  k+i=2,  we  run  into  divergent 
20  11  02 

terms,  namely,  I2Q;  and  although  the  remaining  two  terms  that  do  not 

.  ^  ,  20  ^  02  „  .  ,  , 
vanish,  namely,  and  are  convergent.  Hence,  all  that  we  can  rigorously 

deduce  from  this  expansion  is  the  first-order  correction  to  the  linear  isoprobs. 

From  equations  (4.21),  (4.32),  (4.34),  (4.37)  and  (4.38)  we  are  led  finally  to 


P(A)  « 


^b/Pb 
Pa  Pb 


{ 


^A-P 


2 

A 


P 


A 


J-  2 

U  -Ll 

B  ^B 


B 


(4. 39) 


Of  course,  the  lead  term  is  again  merely  the  familiar  (1-6)  in  a  slightly 
different  form. 

As  a  test  of  (4. 39)  we  use  the  same  example  we  have  considered  hitherto, 
viz.,  (3.27),  (3.29)  and  (3.30);  we  find 

P(A)  of  eqn.  (4.  39)  =  |  +  ^  I  +  2)  =  |  +  ^  =  {|  =  •  72222  .  (4. 40) 
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The  error  here  is^  comparing  this  result  with  (3.30)^  positive  and  precisely 
one  part  in  168^  or^  approximately  0.6%.  This  is^  needless  to  say^  a  great 
improvement  over  any  of  the  approximations  we  have  so  far  derived. 

The  success  of  equation  (4.39)  in  this  case  makes  one  wish  that  one  might 
extend  it  at  least  to  the  next  order,  which  would,  as  we  remarked  below  (4. 38) 
entail  five  new  terms.  The  divergences  already  cited,  however,  seem  for  the 
present  effectively  to  block  this.  Still,  it  is  very  difficult  to  feel  that 
(4.37)  does  not  really  contain  all  the  stochastic  information  inherent  in  the 
duel,  provided  a  suitable  interpretation  of  (4. 36)  is  found.  One  might,  for 
example,  call  upon  Hadamard's  notion  of  the  "finite  part"  of  an  integral;  or 
one  might  attempt  to  confine  attention  to  contributions  of  residues  in  the  right 
half-plane  alone.  At  this  writing,  however,  our  efforts  in  this  direction  have 
come  to  nothing. 

We  close  this  section  with  a  brief  examination  of  a  couple  of  properties 
of  the  of  equation  (4.20).  Substituting  the  latter  into  (4.14)  and  reversing 
the  order  of  summation,  we  find 


v-1 

k=o 


^  [(v-k)}!  ,  -Z]X  ] 

^  v-k-1  v-k 


k 

23k 


k+X 


.k+1 
(l-jiz)  X=o 


V-1 

=  vz  u 


v-1 


(4.  41) 


which  becomes,  when  we  clear  of  fractions. 


E(k)[(v-k)i^_j^_^-zp^_j^]  (l-pz)"'-^-^ 
k=o  X=o 


(4.42) 
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The  highest  power  of  z  which  occurs  in  this  identity  is  the  v — ^  and  when  we 
equate  its  coefficient  on  both  sides  of  the  equation,  we  find 

•  (4.43) 

k=o 

Letting  v=l, 2,3  and  4  in  this,  we  are  easily  led  to  equations  (4.21),  (4.22), 
(4.24)  and  (4.27)  again. 

Consider  next  the  generating  function 

00  \f 

g(x)  =E(-)\^  ^  J  (4.44) 

v=o 


then  (4.43)  yields 


00  v 

EX  V 

— T  l-L  1 

Vl  ^  '^V-1 

V=1 


00  V  V 

M.X  E 


v=o 


v=l  k=o 


v-k^kk 


«  it,,  X  p 

=  Sc-)"  r®  S 


V  v-k-1 


v-k 


^  kl  (v-k)l 

k=o  v=k+l  '  ' 


=  ^  E(-)^  ^  "  ■jg(x)  [cp(px)-l],(4.45) 

^  k=o  v=l  ^ 


from  which  it  follows  that 


"  V  2  ,  ,, 

2i_  =  P, 

^1  cp(px)-l 


v=o 


(4.  46) 


13  September  1963 


29 


SP-1017/003/00 


Here  again  we  readily  check  the  validity  of  the  equations  listed  at  the  end  of 
the  previous  paragraph.  We  note, furthermore,  that  for  an  exponential  firing¬ 
time,  cp(x)  =  l/(l-(j.x)  and  the  right-hand  side  of  (4.46)  becomes  identically 
unity,  confirming  the  remark  made  in  the  paragraph  containing  (4.19). 

We  may  also  examine  similarly  the  behavior  of  the  coefficients  at  the 
"opposite  end,"  so  to  speak.  The  constant  term  in  (4.42)  yields 

-  ''%.l  ' 

k=o 


which  is  trivial,  in  light  of  (4.21).  The  coefficient  of  z  gives,  on  the  other 
hand. 


c=o 

k=o 


(4.  48) 


which  becomes,  by  virtue  of  (4.21), 


or 


v-1  ^  2 

k=o 


v-1 


/V-L 


k  ^^v-k-1  kl 

k=o 


whence^  on  replacing  v  by  v+1^ 


s  V'v-k'^ki 


k=o 


(4.  49) 


(4. 50) 


v+1 


(4.51) 
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Here,  again,  setting  v=0,l,2  and  3,  in  succession,  we  obtain  verification  of 
(4.  21) -(4.  23),  and  (4.25). 

Let  us  now  consider  the  generating  function 


G(x)  =  ^  ; 


v=o 


(4. 52) 


(4. 51)  then  yields 


K  X  X_  /  V+1 

4_/  1, 1  \ 


v-1  «> 


v=o 


vl  v+1  ~  v!  ^  v! 

v=l  v=o 


=  ^  [9(x)-l]  -  ^cp(x) 


00  V  V  00  \|f  CO  M 

V'  2S_  I  Y'  V'  v-k  V 

--  Lj  vl,^^k^^v-k''^kl  kl  ^  (v-k)!  ^ 
ysQ  k*o  k«o  v=k 


00  \lr  ^  00  H 

k  nr  \  f  \ 

=  2j~  X  2-/  —  X  =  G(x)9(x)  , 

k=o  v=o 


(4.  53) 


thus  showing  that 


00  y 


v=o 


1 

X 


_ 1_ 

^  X9(x) 


(4.  54) 


Once  again  we  have  a  verification  of  equations  (4.  21) -(4. 23),  and  (4.25);  and 
this  time  an  exponential  firing-time  causes  the  right-hand  side  of  (4.54)  to 
vanish  identically,  which  is  as  it  should  be. 

There  is  a  last  identity  which  we  shall  point  out.  From  (4.44)  and  (4.46) 


we  have 
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g(x)  2  2  .  . 

X  [1  XCp(M.X) 


(4. 55) 


and  from  (4.52)  and  (4.54)^ 


G(nx)  H  |axcp(iax) 


(4. 55) 


It  follows  that 


-yr  =  1  +  -  G(nx)  , 
g(x)  ^  ’ 


g(x)G((ix)  =  ia[l-g(x)] 


(4.  57) 

(4.58) 


From  this  we  obtain  the  curious  convolution  identity^ 


k 


A=o  ^  ' 

=  0  (k=o) 


k-l. 


(4.  59) 


5.  A  STOCHASTIC  BATTLE 


When^  without  stratagem^ 

But  in  plain  shock  and  even  play  of  battle, 

Was  ever  known  so  great  and  little  loss 
On  one  part  and  on  the  other? 

- King  Henry  V.  (IV,  viii) 


We  assume  that  two  opposing  forces,  A  and  B,  consist  of  a  and  p  men, 
respectively.  All  men  on  side  A  have  the  same  rate-of-f ire,  r^,  and  the  same 
kill  probability,  p^,  and  their  time-to-fire  is  exponentially  distributed. 
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Then  we  have  seen  that  their  time-to-kill  will  be  exponential^  with  rate-of-kill 

equal  to  r  p  .  Making  the  letters  A  and  B  do  double  duty  for  simplicity  of 
£\  £\ 

notation^  we  write 


*  =  Va  ' 


(5.1) 


Since  all  a.  men  on  side  A  are  firing  simultaneously^  the  time  to  A’s  first  kill 
of  one  of  B’s  men  is  the  smallest  in  a  sample  of  size  a  from  an  exponential 
variate  of  mean  1/A.  By  the  well-known  formula  from  order  statistics^  therefore^ 
the  probability  that  A’s  first  kill  occurs  at  time  >  t  is  equal  to 


H^Ct) 


r  -Ati(D:  -Aat 

Le  J  =  e 


(5.2) 


and  so  his  time-to-f irst-kill  is  distributed  as 


u  "Aat 

h^(t)  =  Aae  , 


(5.3) 


i.e.  ^  it  too  is  exponential^  the  rate-of -first-kill  being  merely  A(2;.  A  similar 
result  holds  for  B: 


hg(t)  =  Bpe 


-Bpt 


(5.4) 


Hence^  the  probability  that  the  first  kill  is  inflicted  by  A  on  B  (rather  than 
by  B  on  A)  is^  just  as  in  the  case  of  the  simple  duel  with  exponential  firing- 


times: 
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Prob. 


A  first  kills  a  B  = 


Ag 

AOfBp  ‘ 


(5.5) 


If  we  consider  a  point  in  the  (g, 3) -plane  of  coordinates^  this  means  that 
precisely  two  transitions  are  possible,  and  that  their  respective  probabilities 
are 


Prob{(g,p)  -V  (g,p-l)}  =  ,  (5.6) 

and 

Prob{(a,p)  ->  (g-1,3)}  =  •  (5.7) 

At  this  juncture  it  might  be  argued  that,  once  a  kill  has  occurred  on  one 
side  or  the  other,  the  situation  becomes  very  complicated,  since  we  must  now 
wait  to  see  whether  the  second  kill  for  the  successful  side  comes  before  or  after 
the  still-awaited  first  kill  for  the  other  side  and  the  developing  pattern  of 
events  ramifies  out  more  and  more.  That  this  is  not  the  case  it  is  quite  easy 
to  see.  For,  an  exponential  distribution  describes  a  constant-risk  situation, 
and  is  inherently  Markovian.  Thus,  to  say  that  A’s  time  to  his  first  kill  is 
distributed  as  (5.3)  is  the  same  as  to  say  that  in  every  time  interval  dt  the 
probability  that  A  scores  his  first  kill  is  Agdt,  regardless  of  what  may  have 
happened  prior;  the  only  assumption  necessary  is  that  A  does  have  g  men  during 
the  interval  dt.  Of  course,  if  B  gets  in  first  with  a  kill,  then  A’s  probability 
of  a  kill  immediately  drops  to  A(g-l)dt. 

Thus,  equations  (5.6)  and  (5.7)  fully  characterize  the  course  of  the 
engagement  as  it  unfolds.  They  are  of  course  merely  a  description  of  a  two- 
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dimensional  random  walk  in  which  only  downward  and  leftward  steps  are  permissible 
and  for  which  the  transition  probabilities  are  not  constant  but  vary  from  lattice 
point  to  lattice  point.  Clearly  the  walk  can  never  pass  either  of  the  axes 
since  the  number  of  combatants  on  either  side  is  necessarily  positive.  Now^  in 
practice^  a  military  unit  which  is  suffering  a  sufficiently  high  casualty  rate 
will  break  and  run^  but  for  purposes  of  simplicity  we  shall  assume  that  both 
sides  stick  it  out  until  one  side  or  the  other  is  utterly  annihilated.  It  would 
be  simple  enough  to  treat  the  problem  on  the  assumption  of  a  built-in  breaking 
pointy  but  the  present  discussion  will  be  adequate  for  a  first  treatment.  Thus^ 
the  random  walk  must  with  probability  unity  ultimately  be  absorbed  into  either 
the  a-axis  or  the  p-axis^  the  former  event  representing  a  victory  by  the 

latter  a  victory  by  A.  A  moment's  consideration  will  show  that  a  tie  is 
impossible  since  the  origin  could  only  be  reached  from  either  the  point  (1^0) 
or  the  point  (0^1)^  both  of  which  are  absorbing  states.  We  should  like  to 
determine  the  probability  that  A  wins  and  the  distribution  of  the  number  of  his 
survivors^  given  that  he  wins. 

We  define  p^(o:_,P)  as  the  probability  that  A  wins  starting  with  a  men  opposing 
and  that  he  has  precisely  v  men  left  at  the  termination  of  the  battle.  In 
other  words^ 

p^(a^P)  =  prob.  of  absorption  at  (v^O)^  starting  from  (a^P)  .  (5.8) 

It  then  immediately  follows  from  (5.6)  and  (5.7)  that 
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for  his  path  to  (v^o)  must  pass  through  one  or  the  other  of  the  points  (a-l^P) 
or  (a^p-1).  We  note  also  that 

p^(a,p)  =0  (V  >  a)  ,  (5. 10) 

since  rightward  steps  are  prohibited.  Clearly  the  over-all  probability  that 
A  wins  is  given  by 

a 

P(A)  =Ep^(0!,P)  .  (5.11) 

V  =  1 

We  note  that  here  we  are  using  our  conventional  notation  for  the  simple  duel  in 
using  A  as  an  "argument"  on  the  left-hand  side  of  the  equation;  the  a  and  p  on 
the  other  side  of  the  equation  are  arguments  in  the  customary  mathematical 
usage  of  the  term. 

The  random  variable  v  is  unnormalized^  because  of  the  possibility  of 

absorption  of  the  process  into  the  p-axis.  Even  for  an  unnorraalized  variate^ 

th 

however^  we  can  define  moments  by  means  of  the  usual  formula  for  the  k —  moment 
about  the  origin^ 

=  Z)v*^p^(a,P)  .  (5.12) 

v=l 

The  notation  makes  explicit  the  dependence  of  the  moment  on  the  initial  point. 
Plainly  we  have 
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and  in  order  to  convert  to  moments  in  the  traditional  meaning  we  have  only  to 
divide  (5.12)  by  (5.13).  By  virtue  of  (5.9)  we  now  have 

a  , 

(Ac>^-B^)^l  (a^p)  {Bpp^(a-l,p)  +  Aap^(a,p-1)} 

v=l 

(5.  14) 

,  a-1  ,  a 

=  cTbpp  (a-l,p)  +  bpDv  p^(o!-1,P)  4-  Aa  Sv'^P^(a,p-l)  . 

V=1  V=1 

But  the  first  term  on  the  right  vanishes  because  of  (5.  10)_,  and  (5.12)  evaluates 
what  remains  on  the  right;  we  find 

(AQ!+Bp)pj^(a,P)  =  Bp^j^(a-l,p)  +  AQ!i^(a,p-l)  .  (5.15) 

This  equation  is  of  precisely  the  same  form  as  (5.9)^  with  the  subscript  v 

dropped.  It  is  interesting  to  remark  that  (5. 15)  does  not  depend  at  all  upon 

k  itself.  This  is  not  to  say,  of  course,  that  M-j^(0!,  P)  is  independent  of  k, 

since  the  latter  enters  into  the  boundary  conditions. 

In  this  respect  we  note  that  specifying  the  values  of  along  the  (positive) 

a  and  p  axes  serves  to  determine  the  solution  of  (5. 15)  uniquely  throughout  (all 

the  lattice  points  of)  the  first  quadrant.  For,  all  we  need  to  know,  by  (5.15), 

to  find  the  value  of  at  a  given  point  are  its  values  at  the  lattice  points 

k 

next  left  and  next  below.  We  can  therefore  build  up  the  solution  from  left  to 
right  a  tier  at  a  time.  As  for  the  appropriate  boundary  conditions  we  note 
that  (5.8)  patently  implies  that 


p^(a,0)  =  1  (v=0!)  , 

=  0  , 


(5.  16) 
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and 


P/0,P)  =  0  (all  V)  .  (5.17) 

Putting  these  values  into  (5.12)^  we  obtain 

Pj^(a,0)  =  ,  (5.  18) 

and 

Pj^(0,p)  =  0  .  (5.  19) 

In  a  recent  paper,  Brown  has  considered  the  special  case  where  k=0;  in 
other  words  he  has  confined  his  attention  to  consideration  of  P(A).  He  obtained 
for  it  an  expression  which  we  unknowingly  duplicated;  the  reader  is  referred  to 
his  paper  in  the  bibliography  for  the  details.  The  expression  follows  readily 
enough  from  successive  applications  of  (5.15),  together  with  (5.18)  and  (5.19); 
it  is  an  alternating  series  containing  exactly  a  terms  and  is  of  limited 
practical  applicability  for  large  a  because  the  individual  terms  become  enormous 
and  a  very  great  cancellation  of  significant  figures  results,  so  that  in  order 
to  obtain  meaningful  answers  one  would  have  to  carry  an  inordinate  number  of 
decimal  places.  Equation  (5.15)  itself  involves  no  such  loss  of  accuracy  and 
is  quite  suitable  for  a  and  p  not  too  large,  and  could  very  readily  be 
programmed  for  a  digital  computer.  Here  again,  however,  for  increasing  a  and 
P  one  finds  that  the  computation  time  mounts  rapidly.  What  we  require  is  a 
serviceable  asymptotic  expansion.  We  shall  now  present  merely  the  zeroth-order 
behavior  of  the  solution  for  large  OL  and  P;  we  have  not  yet  succeeded  in 
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extending  the  analysis  any  further  (notwithstanding  repeated  attempts)^  but 
the  result  which  emerges  seems  to  us  interesting^  unexpected^  and  indicative  of 
the  desirability  of  pursuing  the  matter. 

We  note  first  the  Taylor’s  expansions 

+  ~  ...  y  (5.20) 

and 

+  “  •••  •  (5.21) 


Since  we  are  looking  for  a  "stabilized"  asymptotic  behavior  for  large  a  and 
it  is  reasonable  to  ignore  the  higher-order  terms  here;  the  two  Taylor’s  series 
are  in  powers  of  unity^  which  we  assume  small  compared  with  a  and  p.  Substituting 
into  (5.  15)^  we  find 


(5.22) 


we  here  suppress  the  arguments  (CC^P)  since  they  are  the  same  as  the  independent 
variables.  Now^  the  subsidiary  equations  necessary  for  the  solution  of  (5.22) 
--  cf.  Forsyth  --  are 


da  _  dB  ‘^^k 
Bp  Aa  “  0  ^ 


(5.23) 


the  last  two  of  which  yield 
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being  the  constant  of  integration.  Similarly^  the  first  two  equations  yield 


Aada  =  Epd^  ,  Aa^-Bp^  = 


(5.25) 


Hence  the  most  general  solution  of  (5.22)  is 


C,  =  ^  (CJ  , 


(5.26) 


where  ^ is  a  completely  arbitrary  function,  soon  to  be  specified  by  the  boundary 
conditions.  Employing  (5.24)  and  (5.25),  (5.26)  becomes 


(5.27) 


But,  from  (5. 18)  we  have 


'if  (Aa^)  =  . 


(5.28) 


When  we  write 


u  =  Aa^  ,  a  =  (u/A)^^^  , 


(5.29) 


(5. 28)  becomes 


'if  (u)  =  (u/A)’"/2  ^ 


(5.  30) 


and  using  this  in  (5.27)  we  have,  finally, 


,  ,2  B  2.k/2 

(I.  (a,p)  =  (a  -  7  p  ) 


(5.31) 
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These  are^  of  course,  moments  about  the  origin,  and  the  cases  k=l  and  k=2 
provide  us  with  the  mean  and  variance,  namely, 

H(a,p)  =  (5.32) 

and 


a^(a,p)  =  =  0  . 


(5.  33) 


It  is  the  last  result  which  comes  as  a  surprise,  but  looking  at  (5.31)  we  see 
that  it  is  indeed  precisely  the  situation  which  obtains  when  we  have  a  "random" 
variable  which  is  in  fact  constant  and  therefore  identically  equal  to  the  value 
given  by  (5.32).  Now,  it  is  clear  that  for  (5.32)  to  be  a  meaningful  statement 
we  next  have 


When  this  holds,  then  (5. 13)  and  (5. 31)  show  that 


(5.  34) 


P(A)  =  1  .  (5.35) 

In  other  words  in  the  wedge  between  the  a  axis  and  the  line  0!  =  X  ^  ^  ^  always 
wins  (assuming  that  a  and  p  are  both  large)  and  he  always  finds  himself  at  the 
end  of  the  fray  with  precisely  the  number  of  survivors  given  by  (5.32).  Of 
course,  a  complementary  situation  obtains  in  the  remainder  of  the  first  quadrant. 

We  remark  that  (5.32)  is  exactly  what  we  get  from  Lanchester's  "square" 
law  (Morse  and  Kimball)  if  we  set  A's  efficiency  relative  to  B  equal  to 
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E  =  A/B  (5.36) 

and  put  the  number  of  combatants  on  side  B  equal  to  zero.  Thus^  we  have 
established  that  Lanchester ’ s  square  law^  which  is  based  on  strictly  deterministic 
reasoning^  is  in  fact  valid  and^  furthermore^  valid  in  a  strictly  deterministic 
sense  (i. e. ,  with  no  random  component  to  a  first  approximation)^  when  large 
numbers  face  one  another.  But^  more  than  that  we  have^  by  (5.1)  and  (5.36)^ 
been  able  to  identify  his  relative  efficiency  in  a  very  simple  and  natural 
manner. 

It  is  of  interest  to  remark  that  the  solution  of  (5.22)  calls  upon  only 
one  of  the  boundary  conditions,  namely,  (5.18),  whereas  the  solution  of  (5.15) 
calls  upon  both  (5.18)  and  (5.19).  It  is  this  that  accounts  for  the  discontinuity 
in  the  solution  of  the  former.  Equation  (5.35)  is  less  anomalous  when  we 
consider  that  (5. 18)  implies  that  P=1  along  the  a-axis,  and  hence  bV/bk  =  0 
there.  But,  setting  [3=0  in  (5.22)  shows  that  oP/op  =  0  also,  along  the  a-axis. 
Thus  P  must  be  constant  and  hence  equal  to  its  value,  unity,  along  the  a-axis. 

The  overwlielming  effect  of  superiority  of  numbers  has  long  been  recognized,  and 
reference  to  the  index  of  Clausewitz's  On  War  will  show  the  reader  that  that 
work  is  laced  through  with  injunctions  as  to  the  importance  of  this  factor. 

The  author  inveighs  against  quixotic  generals  who  have  deliberately  failed  to 
bring  the  full  force  of  their  troops  to  bear  upon  the  enemy,  although  as  a 
practical  man  of  affairs  he  occasionally  tempers  his  position  with  a  warning 
against  the  over-simplification  of  assuming  that  there  are  no  other  factors 
involved  in  the  outcome  of  a  battle  besides  superiority  in  numbers. 


The  fact  that  k  does  not  appear  explicitly  in  (5. 15)  means  we  can  rewrite  that 
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equation  very  simply  in  terms  of  tli  tte  moment-generating  function.  Define 

“  k  )|( 

cpCzjO^P)  =  u'T7'  (5.37) 

k=o  " 

Then  (5. 15)  yields  at  once 

(Aa+Bp)cp(z;Q;,p)  =  Bpcp(2;i;s2;(l-  1,  p)  +  AQtp(z;ci;,p-l)  ,  (5.38) 

and  the  approximate  analysis  folloilloss  through  exactly  as  before^  leading  to  the 
counter-part  of  (5.27),  viz., 

cp(z;a,p)  =  'i^(Aa^-BpV.  .  (5-39) 

The  boundary  condition  is  now,  by  [oyvaLr  tue  of  (5.18)  and  (5.37),  given  by 

cp(z;a,0)  =  Zy  ^  *  n  .  (5.40) 

k=o 


In  equation  (5.38),  z  plays  the  rS  iSle  of  a  parameter,  just  as  k  did  in  (5.15). 
Equations  (5.39)  and  (5.40),  togetJigetlier  with  (5.29)  now  yield 


cp(zja,p) 


(5.41) 


thereby  checking  (5.31).  It  seetiasis  to  be  pretty  much  of  a  toss-up  whether  this 
analysis  or  the  former  is  any  siji|iiiiipler  in  this  case,  although  this  does  not 
preclude  the  possibility  that  oneii«e(>:f  them  might,  in  a  more  fully-blown  treatment, 
turn  out  to  be  preferable  to  the  oae otlner. 

We  assumed  at  the  beginning  ogg of  Chis  section  that  firing-times  were 
exponential  and  noted  in  (5.  3)  thiiJttat  this  implied  that  time-to-f  irst-kill  was 
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also  exponential.  It  makes  sense  to  ask  what  happens  when  we  drop  this  restriction. 
The  extreme-value  distribution,  of  the  algebraically  smallest  value  in  a  sample 
of  large  size,  is  known  to  follow  a  Fisher-Tippett  distribution  (cf.  Cramer), 
provided  certain  conditions  are  met,  one  of  which  is  that  the  distribution 
being  sampled  from  have  an  infinite  tail  in  the  negative  direction.  This  is 
plainly  not  the  case  with  the  inherently  positive  firing-times  we  have  to  deal 
with,  and  we  may  expect  that  result  to  break  down;  it  does.  For  concreteness 
let  us  take  the  time-to-fire  to  be  distributed  as  a  gamma-variate  with  parameters 
2  and  v; 


f(t) 


(5.  42) 


This  will  obviously  be  a  more  realistic  assumption  than  the  exponential.  It 
follows  that 


F(t)  = J  f(t)dt  =  (l+vt)e 


•vt 


(5.43) 


upon  integrating  by  parts.  Hence  the  right-sided  cumulative  distribution  function 
of  the  time-to-f irst-kill  is,  analogously  to  (5.2), 


^  /■,  nO  -ovt 

H(t)  =  (1+vt)  e  , 


(5.  44) 


where  a  is  still  the  number  of  men  on  side  A.  Thus,  the  frequency  function  is 


h(t)  =  -H>(t)  . 


(5.45) 
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Now^  for  any  positive  random  variable^  and,  in  particular,  in  the  present  case, 
the  mean  is  given  by 

00  00  00  00 

p  =  ThCOtdt  =  -  /  H>(t)tdt  =  -  tH(t)  +  /  H(t)dt 
j  o  O  -^0^0 

00 

=  rH(t)dt  ,  (5.46) 

J  o 

provided  H(t)  goes  to  zero  as  t  ->•  “  faster  than  1/t,  i.  e. ,  provided  that 


H(t)  =  o  (1/t)  (t  -»-  “)  . 


(5.47) 


Therefore  the  mean  time-to-f irst-kill  is,  in  the  case  sub  judice , 


p  =  /  (1+vt)  e  dt  , 


(5.  48) 


on  account  of  (5.44)  and  (5.46).  Consequently, 


o  k=o 


00 

k.k  “Qyt  V'/Ct  k  r  Jf.  -Qytj 
2J(pv  t  e  dt  =Z7(pv  /  t  e  dt 


.  V  A  ^  kl  _  ^  ^  1 

Zy^k^^  k+1  ocv  ^  ,  i\i  ^ 

k=o  (av)  k=o  (a-k)  la 


,  a  k 
al 

ofl  kl  ^ 
va  k=o 


(5.  49) 


upon  reversing  the  order  of  summation.  Now,  the  terms  in  the  summation  are 
proportional  to  the  probabilities  in  a  Poisson  distribution  of  parameter  a. 

This  has  mean  a  and  standard  deviation  and  approaches  normality  ’when  a  ■>  “. 
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The  terms  for  k=a-l  and  k=0!  are  readily  seen  to  be  equal^  and  the  last  argument 
shows  that  the  totality  of  terms  contributing  appreciably  to  the  sum  lie  within 
a  few  multiples  of  •Ja  on  either  side  of  these  terms.  Thus^ 


“  c'' 


Z\ 

'  kl 

K=0 


k 

kl 


00 

^k=o 


1  a 


(a  ->■ 


(5.  50) 


Substituting  this,  and  Stirling’s  formula,  into  (5.49)  we  get 


jj.  ^ 


1  -a-1  rr-  -ala,  . 

—  a  -a  ‘=^v2jre  ‘  2  ^  '*  “)  > 


(5.51) 


which  is  to  say  that 


^  /  2a  ■  V  <"  ^  ”)  • 


(5.52) 


We  therefore  introduce  a  new  random  variable. 


T  2  'Ja  vt  , 


(5.53) 


for  which  the  expectation  is 


E(t) 


^/T 


(a  ->  «>) 


(5.54) 


Its  RSCDF  (right-sided  cumulative  distribution  function)  is,  by  (5,44), 


H(V  =  (1  . 


(5.55) 


which  yields 
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log  H  =  -  \^oa  +  a  log  (1  + 


/  /  1  a.  ^  4.  ^ 


(a  -»■  “) 


(5.56) 


Thus^ 

H(T)  .V.  exp(-  i  T^)  (Q!^“)  ,  (5.57) 

and^  by  (5.45)^ 

1  2 
■  2 

h(T)  -V  Te  (a  -»■  “)  .  (5.58) 

It  is  a  simple  matter  now  to  check  the  validity  of  (5.54). 

We  see  from  this  last  analysis  that  the  time-to-f irst-kill  behaves  quite 
differently  from  the  way  it  did  before;  contrast  (5.58)  and  (5.53)  with  (5.3). 
However^  although  the  battle  commences  differently  in  this  case^  we  may  legiti¬ 
mately  anticipate  that  it  will  soon  shift  over  to  the  situation  described  earlier. 

For^  consider  first  of  all  a  men  firing  at  a  target^  each  with  mean  time-to- 

2 

fire  p  and  variance  a  .  Then  the  a  first  rounds  fired  by  them  will  produce  a 

near-volley  which  is  a  statistical  image  of  the  firing-time;  the  a  second 

rounds  fired  will  be  spaced  out  a  little  more  perceptibly^  having  mean  2p  and 

variance  2a^;  and  so  on.  The  n—  round  will  thus  be  a  sample  of  size  a  from 

a  distribution  with  mean  np  and  variance  no"^;  and  the  (n+1)—  round  will  have 

2 

mean  (n+l)p  and  variance  (n+1) a  .  Thus  successive  rounds  will  begin  to  overlap 
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and^  provided  a  is  not  too  small^  this  overlapping  will  ultimately  result  in  a 
complete  jumbling  of  the  rounds^  so  that  the  timing  of  the  successive  shots  as 
they  actually  occur  will  degenerate  into  a  Poisson  process.  The  fact  that  B  is 
simultaneously  firing  at  A  is^  to  be  sure^  an  added  complication^  but  it  does 
not  alter  the  conclusion  that  (5.6)  and  (5.7)^  with  A  and  B  suitably  intepreted^ 
will  be  valid  for  the  bulk  of  the  process. 


6.  ANOTHER  STOCHASTIC  BATTLE 

Or  what  king^  going  to  make  war  against  another  king^  sitteth 
not  down  first^  and  consulteth  whether  he  be  able  with  ten 
thousand  to  meet  him  that  cometh  against  him  with  twenty 
thousand? 

— Luke  14:  31 

It  is  evident  that  the  difficulties  intrinsic  to  the  discussion  presented 
in  the  last  section  are  immediately  traceable  to  the  fact  that  the  transition 
probabilities  (5.6)  and  (5.7)  are  functions  of  a  and  p.  Let  us  therefore 
replace  these  equations  by 

Prob[(a^P)  (a^p-1)}  =  p_,  (6.1) 

and 

Prob[(a^P)  ^  (a-l,P)}  =  q  ,  (6.2) 


where  p  and  q  are  constant  and^  of  course^ 
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p  +  q  =  1  .  (6.  3) 

This  battle  may  be  conceived  of  as  a  series  of  single  contests  between  A's  and 
B’s  men^  there  being  constant  probability  p  that  A  wins  any  one  of  these  contests. 
A  graphic  way  of  bearing  in  mind  the  difference  between  the  two  battle  schemata 
is  as  follows:  in  the  former  battle  the  opposing  forces  were  deployed  facing 
one  another  along  parallel  lines  and  all  surviving  men  on  both  sides  were 
engaged  in  conflict  at  all  times;  in  the  present  battle  the  opposing  forces  are 
deployed  in  two  columns  which  meet  head-on  in  a  narrow  defile  and  combat  is 
confined  to  the  men  heading  the  two  columns  of  survivors.  As  before  we  were 
led  to  a  justification  of  Lanchester’s  square  law^  here  we  might  well  expect  to 
vindicate  his  linear  law,  and  this  we  shall  indeed  find  to  be  the  case. 

We  retain  the  notation  of  the  last  section  and  begin  by  skimming  through 
those  parts  of  the  analysis  which  carry  over  with  trivial  alterations.  Equation 
(5.9)  becomes 

p^(a,p)  =  qp^(a-l,p)  +  pp^(a,p-l)  ,  (6.4) 

which  leads  to 

+  P^^j^(o:,P-l)  ,  (6.5) 

analogous  to  (5. 15),  and 

q  ^  +  P  ^  =0  (a  and  p  both  large)  ,  (6.  6) 

analogous  to  (5.22);  all  the  rest  of  the  intervening  equations  (except,  of 
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course;  for  (5.  14))  remain  unchanged. 


q  p  0  ^ 


The  subsidiary  equations  (5.23)  now  become 

(6.  7) 


yielding  (5.24)  and;  in  place  of  (5.25); 


pa  -  qp  =  . 

Hence  the  most  general  solution  of  (6.  6)  is 


(pa  -  qp)  ; 

analogous  to  (5.27).  Boundary  condition  (5.18)  now  yields 


^  (pa)  =  a*^  ; 

or;  setting 


u=pa;  a=u/p; 
^  (u)  =  (u/p)^  . 


(6.8) 


(6.9) 


(6.  10) 


(6.  11) 
(6.  12) 


Substituting  these  results  into  (6.9);  we  find 

Hl,(a;P)  =  (a  -  ^  p)^  ;  (6.13) 

which  implies  that  the  random  variable  v  implied  in  (6. 4)  is  in  fact  identically 
equal  to  a  -  P 
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V  H  a  -  P  (a  and  p  both  large) 


(6. 14) 


Thus^  analogous  to  (5.34)  and  (5.35),  we  have 


P(A)  =  1 


(a  >  -^  p)  . 


(6. 15) 


Hence  A  always  wins  over  B  provided  a  >  -^  P  (assuming  a  and  p  both  large)  and 
his  victory  costs  him  precisely  p  of  his  men.  If  we  set  the  number  of  men 
on  side  B  equal  to  zero,  and  here  define  the  exchange  rate  to  be 


E  =  p/q 


(6. 16) 


we  are  led  to  (6.14). 

But,  we  said  at  the  outset  that  we  might  expect  a  detailed  analysis  to  be 

more  tractable  for  this  schema  than  for  the  former  and  we  shall  now  see  that  this 

is  perfectly  true.  For,  the  process  may  now  be  visualized  as  a  two-dimensional 

random  walk  with  downward  transition  probabilities  all  equal  to  p  and  leftward 

transition  probabilities  all  equal  to  q.  For  A  to  wipe  out  B,  with  a  loss  of 

precisely  v  of  his  men,  requires,  when  he  starts  from  the  point  (a, P),  that  he  take 

exactly  v  steps  to  the  left  and  exactly  p  steps  down,  which  introduces  a  factor 
V  6 

of  q  p^.  (To  be  sure,  the  order  in  which  these  leftward  and  downward  steps  are 
taken  will  vary  from  one  realization  to  the  next,  but  the  total  numbers  of  the 
two  different  kinds  of  steps  clearly  must  always  be  v  and  p,  respectively.) 

Now,  of  these  steps,  the  nature  of  the  last  alone  is  dictated;  it  must  be 
downward,  from  (v,l)  to  (v,0),  for,  since  the  a-axis  is  an  absorbing  barrier,  the 
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step  (v+1^0)  to  (v^O)  is  prohibited.  Hence  the  total  number  of  permissible  paths 
is  the  number  of  ways  that  v  leftward  steps  can  be  chosen  out  of  a  total  of 
P+v-1  steps.  We  therefore  see  that 


>P+v-l.  V  p 
=  (  V  )9  p 


(6.  17) 


It  is  a  simple  matter  to  check  that  this  satisfies  the  recurrence  (6.4)  and  the 
boundary  conditions  (5.16)  and  (5.17).  It  follows  from  (6.17)  that 


P(A) 


a-1 

S( 

v=o 


p+v 

V 


-1.  V 

)q 


(6.18) 


and  combining  this  with  the  corresponding  result  for  P(B)  we  are  led  to  the 
interesting  and  non-trivial  identity 

p*"  +  <>“  ‘"srrbp"  - 1  .  (6. 19) 

v=0  v=0 


For  the  purpose  of  clarity^  solution  (6.  18)  is  set  forth  explicitly  for  the 
first  few  cases  in  the  following  table. 


3 

0 

3 

P 

p^(l+3q) 

(l+3q+6n^) 

2 

0 

2 

P 

p^(l+2q) 

p‘'(l+2q+3q") 

1 

0 

P 

p(i+q) 

p(l+q+q  ) 

0  J 

lE 

1 

1 

1 

cr-^O 

1 

2 

3 

Table  2.  P(A)  qua  function  of  a  and  p. 
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The  symbolic  argument  A  in  P(A)  is  of  no  utility  analytically  and  so  henceforth 
we  write 


P(a,p)  =  P(A)  , 


(6.20) 


understanding  that  the  probability  of  win^  P.  refers  to  side  A,  and  displaying 
explicitly  its  dependence  upon  the  initial  numbers  of  forces.  Reading  up  along 

any  given  column  in  Table  2  we  see  that_,  for  fixed  aj  P(q:,3)  goes  from  1  to  0  as 

P  goes  from  0  to  «>;  in  other  words  it  behaves  like  a  RSCDF  with  p  as  RV  (random 
variable).  This  suggests  that  we  might  hope  to  get  a  firmer  hold  on  (6.15)  by 

looking  at  the  moments  of  P;  thus^  we  would  anticipate  that  its  mean  should  fall 

somewhere  near  pa/q^  and  that,  since  P  apparently  drops  precipitately  thereabouts, 
the  standard  deviation  of  p  ought  to  be  small,  relative  to  its  mean.  We  proceed 
to  follow  up  this  clue. 

For  a  discrete  positive  RV,  p,  with  frequency  distribution,  f^,  we  defii 
the  RSCDF  to  be 


line 


Fp  ^  Sfb 
P  b=p  ^ 


(P=0,l,2,...)  , 


(6.21) 


which  makes 


F  =  1 
o 


(6.22) 


and 


f  =  F  -  F 

P  P  P+1 


(P— 0 ,1,2,...)  . 


(6.23) 
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Thus 


the  k- 


th 


factorial  moment 


is 


E[p(p-1)-  •  •  (p-k+1)]  =  Sp(e-l)-*  •  0-k+l)f 

3=0 

00 

=  Ep(P-I)  -  •• 

fi=k 


00  00 
=  Sp...  0-k+l)F  -  E  (B-1)---  (p-k)F 
;3=k  ^  p=k-l  l  ' 


=  Eo-1)---  (e-k+i)tp-o-k)]F 

p=k 


=  k  E(P-l)(P-2)---0-k+l)F„  (P=l,2,3,...)  .  (6.24) 

p^-k  ^ 


We  havc^  in  the  discussion  preceding  (5.8)^  remarked  the  inaccessibility 
of  the  point  (0,0)  to  the  recurrence  equation^  a  fact  which  has  already  been 
duly  indicated  in  Table  2.  We  now  want  to  set  up  a  generating  function  for 
P(a^P),  and  so^  keeping  this  difficulty  rn  mind^  we  write 


g(>^^y)  =  Z)y^  ^p(a^e)  •  (6.25) 

0=1  p=l 


The  reason  for  the  exponent  on  y  is  that  it  makes 


0=1 


Z;(P-l)(P-2).-.  (p-k+l)P(o,p)  ; 
P=1 


(6.26) 


the  first  k  terms  in  the  inner  summation  vanish  and,  in  consequence,  (6.  24)  and 
(6.26)  yield 
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E[pO-l)-.  .  O-k+1)] 


where  the  notation  introduced 


if  f(x)  =  IjC  X  ^  ^ 
Q?»0  “ 


■ ''  “S"  ,  ’ 

X  Jy=l 

carries  the  construction: 

then  f(x)  =  . 


(6.27) 


(6.  28) 


It  follows  from  (6.18)  and  (6.25)  that 

.  .  a  •A  p-1  p  ^Y^.p+v-1  V 

g(x,y)  =  Zvx  Ijr  ZvC  V  ‘''i 

Q?=l  p=l  V=0 

-  EpV-‘ !;x“  s>r .  (^-25) 

P=1  0=0  v=o 


since^  P  being  always  greater  than  0  now^  the  inclusion  of  the  term  for  Q^O 
introduces  only  a  zero  into  the  summation  over  a.  Inverting  the  order  of  summation 
in  the  innermost  pair  of  sums^  (6.29)  becomes 

g<x,y)  =  Bx“ 

P=1  v=o  0?=v+l 


X  p  p-1  V'/.P+V-ls  ^  ^ 

y  Z/C  V  )9  ^  ; 

P=1  v=o 

which  reference  to  (3. 10)  enables  us  to  write  as 


8(^;y)  =  1?;^  Z/p^y^'^  Z/(-)’^(  v^)q''x’^ 

P=1  v=o 


(6.  30) 
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X 

1-x 


1 

(l-qx)P 


X  p  1 

1-x  1-qx  ,  py  ^ 
1-qx 


(6.31) 


or, 


g(x,y) 


_ E2 _ 

(l-x)[l-(qx+py)] 


(6. 32) 


4 


It  is  now  easy  to  see  that 


=  £X_  (k-1)  Ip*"'^ 

y=l  [l-(qx+py)  Jy=l 


(q-qx)^  (l.x)’^^' 


(6.  33) 


and^  expanding  this  into  a  binomial  series^  we  obtain 


pNk-l  -| 

o  _£ 


dy 


k-1 


y=l 


=  (k-l)i(|)^  hi-f 


Of=o 


=  (k-l)l(^)  E 
^  0=0 


(6.  34) 


upon  yet  another  invocation  of  (3.10).  Equations  (^6.27)  and  (6.34)  now  yield 


E[p(p-i)...(p-k+i)]  =  kKl)*" 


=  k !(£)''  , 


(6.  35) 
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when  we  call  to  mind  the  fact  that 


(v)  =  (n-v)  . 


(6.36) 


Therefore  we  have,  finally, 


E[pO-l)--- O-k+1)]  =  aCcH-l)...  (CH-k-lX^) 


(6.  37) 


The  first  significant  case  here  is  that  for  which  k=l;  we  find 


EO)  =  -^  a  . 

q 


(6.  38) 


This  is  really  rather  more  than  we  dared  desire,  since  it  says  that  the  handicap 
point  described  by  (6.15),  which  marks  a  cross-over  from  a  victory  by  A  to  a 
victory  by  B,  and  which  was  derived  on  the  assumption  that  a  and  p  are  both  large 
holds,  in  a  mean  sense,  for  all  a.  and  p,and,  furthermore,  holds  precisely. 

Letting  k=2,  we  see  from  (6. 37)  that 


E(p^-p)  =  (f)  a(a+l)  , 


(6.  39) 


whence 


E(p^)  =  (^)  a(a+l)  +  ^  a 


(6.  40) 


Thus,  the  variance  of  p  is  given  by 


var  (p)  =  (^)  a(a+l)  +  ^  a 


(I)  a 


2  2 
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Equations  (6.38)  and  (6.41)  are  very  informative,  indeed,  since  they  begin  to 
suggest  that  P  may  actually  be  the  sum  of  Q!  independent  selections  of  the  same 
random  variable.  Upon  this  hint  we  speak;  if  this  is  truly  the  case,  then  the  RV 
being  sampled  can  only  be  that  for  which  Q?=l,  i.  e.  ,  that  whose  RSCDF  is  the  first 
column  of  Table  2.  But  that  is  of  course  merely  a  RV  distributed  according  to  a 
geometric  distribution  whose  frequency  distribution  is,  by  (6.23)  and  (6.18)  with 
Q^l, 

fp  =  p^  -  p^"^^  =  q-p^  .  (6.42) 

Now,  the  factorial-moment-generating  function  of  any  RV  p  is  by  definition 

;t(z)  =  E(z^)  ,  (6.43) 

which  reduces,  in  this  case,  to 

00 

«o(z)  =  q  ZyP^z^  =  •  (6.44) 

P=o  ^ 

But  plainly  the  MGF  with  argument  log  z  is  in  general  given  by 

cp(log  z)  =  E(e^  ^)  =  E(z^)  =  rt(z)  ,  (6.45) 

according  to  (6.43).  The  well-known  properties  of  the  MGF  show  that  if  p  is 
the  sum  of  OC  independent  RV’s,  each  of  which  has  MGF  cp^(z),  i:hen  the  MGF  of  p 
is 


cp(z)  =  [cp^(z)]°^  . 


(6.  46) 
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Calling  the  FMGF  (factorial-moment-generating  function)  of  p  Jt(z),  we  thus  have 

rt(z)  =  cp(log  z)  =  [cp^(log  z)]*^  =  [rt^(z)]“  j  (6.47) 

i. e. ^  the  FMGF  has  exactly  the  same  properties  as  the  MGF.  Here  rt^Cz)  is  the 
FMGF  of  the  RV  being  sampled^  and  the  hypothesis  we  are  entertaining  in  the 
present  paragraph  assumes  that  this  is  given  by  (6.44).  Hence  the  conjecture 
before  us  states  in  effect^  by  virtue  of  (6.44)  and  (6.47),  that  the  FMGF  of  p 
is 


(6.  48) 


Now,  it  follows  from  (6.43)  that 


Ic  1 

^  I  H  :t^^^(l)  -  E[p(p-1)...  (p-k+1)] 
dz  Jz=l 


(6.  49) 


Applying  this  to  (6.48)  produces  the  result 


(1-pz) 


k  -1 


Jz=l 


a  k  k 

=  «(«+!)  •  •  •  (a+k-1)  — =  a(c«-l)  •  •  •  (QH-k-l)  (■^)  ,  (6.  50) 

(l-p)®^*"  ^ 

which  is  identical  with  (6.37).  The  conjecture  is  therefore  established;  namely, 
to  state  it  explicitly:  P(a, p)  is  the  RSCDF  of  a  RV  which  is  the  sum  of  a 
independent  selections  of  a  variate  geometrically  distributed  according  to  (6.42). 
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This  is  truly  rather  a  remarkable  result^  since  it  affirms  that  the  device 
introduced  just  following  (6. 20)  seemingly  as  a  deus  ex  machine  actually  leads 
to  an  exact  and  rigorous  conclusion.  By  (6.45)  and  (6.48),  the  MGF  of  ^  is 
equal  to 


cp(2)  =  rt(e^)  =  (—3—)“  , 


(6.51) 


1-pe 


and  so  its  CGF  (cumulant-generat ing-function)  is,  by  (2.3), 


6(2)  =  a  log  ( — ^)  =  a| 
1-pe 


-z+log  (-^— ) 
1-pe  J 


(6.52) 


from  which  we  can  identify  the  successive 
of  (2. 3).  But,  we  have  come  full  circle, 
necessary  when  we  now  call  upon  (2.5)  and 
to  interchange  p  and  q  and  to  write  z  for 


cumulants  from  the  right-hand  side 
for  no  calculations  whatsoever  are 
(2.22)1  To  do  so  it  is  only  necessary 
6;  we  find,  from  (6.52), 


K,  =  <(-l  +  -)  =  a 
1  q  q 


(6.  53) 


K 


2 


(6.54) 


and 


^3  "  q^^^ 

q 


(6.55) 


the  first  two  of  these  equations  agreeing  with  (6.38)  and  (6.41).  The  last 


two  give  the  skewness. 
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or, 


(6.56) 


(6.57) 


The  proposition  immediately  following  (6.50),  together  with  the  central-limit 
theorem^  enables  us  to  state  that  the  variate 

1  =  ^^  'vN(0,l)  (a^~)  ,  (6.58) 

i.  e.  ,  it  is  distributed  as  a  standardized  normal  variate  for  large  a.  Since  P 
is  the  RSCDF  of  p,  the  last  result  shows  that  if  p  exceeds  pa/q  by  appreciably 

more  than  Va,  then  A  is  virtually  certain  to  lose,  whereas  if  p  falls  short  of 

pa/q  by  appreciably  more  than  n/q:,  then  A  is  virtually  certain  to  win.  The 
tangible  values  may  be  read  from  the  tables  of  the  normal  distribution. 

These  arguments  serve  to  validate  and  delimit  (6.  15)  and^  in  so  doing,  to 
add  increased  credibility  to  (5.34)  and  (5.35),  which  were  derived  by  a  completely 
analogous  chain  of  reasoning.  The  way  is  now  open  for  a  thorough -going 

examination  of  the  "linear"  battle,  along  the  lines  just  initiated;  we  hope  to 

be  able  to  present  further  results  in  a  sequel. 


7.  SUMMARY 


We  have  examined  in  this  report  a  number  of  aspects  which  go  to  make  up 
a  stochastic  theory  of  duels.  We  have  considered  the  manner  in  which  a 
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marksman’s  time  to  hit  a  passive  target  may  be  statistically  characterized  in 
terms  of  his  time  to  fire  a  single  round^  and  have  used  this  information  to 
obtain  approximate  solutions  to  the  duel.  We  have  worked  out  one  case  of  some 
generality  in  which  the  solution  of  the  duel  may  be  written  down  in  a  simple 
form.  We  have  described  two  battle  schemata  which  lead  in  the  limit  of  large 
numbers  of  opposing  forces  to  Lanchester’s  square  and  linear  laws;  for  these 
we  have  broadened  the  stochastic  output  to  include  not  only  the  probability  of 
a  given  side’s  winning,  but  also  the  cost  of  victory  in  the  number  of  his 
forces  sacrificed.  We  have  found  the  latter  schema  to  be  the  mathematically 
more  tractable  situation  and  have  found  there  a  convenient  and  rather  remarkable 
expression  for  the  probability  that  a  given  side  wins,  noting  that  the  techniques 
there  employed  give  promise  of  providing  a  full  solution  to  that  duel. 
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